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In the full nonlinear cosmological perturbation theory in the leading order of the gradi- 
ent expansion, all the types of the gauge invariant perturbation variables are defined. The 
r^ I metric junction conditions across the spacelike transition hypersurface are formulated in a 

manifestly gauge invariant manner. It is manifestly shown that all the physical laws such 
O ' as the evolution equations, the constraint equations, and the junction conditions can be 

\^ • written using the gauge invariant variables which we defined only. Based on the existence 

^ ! of the universal adiabatic growing mode in the nonlinear perturbation theory and the p 

'p'l philosophy where the physical evolution are described using the energy density p as the 

O I evolution parameter, we give the definitions of the adiabatic perturbation variable and the 

'^ • entropic perturbation variables in the full nonlinear perturbation theory. In order to give 

^ c^^ ■ the analytic order estimate of the nonlinear parameter /^vl , we present the exponent eval- 

uation method. As the models where f^i changes continuously and becomes large, using 
the p philosophy, we investigate the non-Gaussianity induced by the entropic perturbation 
of the component which does not govern the cosmic energy density, and we show that in 
order to obtain the significant non-Gaussianity it is necessary that the scalar field which 
supports the entropic perturbation is extremely small compared with the scalar field which 
[-^ . supports the adiabatic perturbation. 
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§1 Introduction and summary 

In the inflationary scenario, the quantum fluctuations of the scalar fields driving the in- 
flationary expansion of the universe are the origins, that is the seed perturbations of the 
temperature fluctuations of the cosmic microwave background radiation (CMB) and the 
cosmic large scale structures such as galaxies and clusters of galaxies. These seed pertur- 
bations generated in the horizon during the inflationary expansion are stretched and go 
out of the horizon. They stay outside the horizon until they return into the horizon in the 
Friedman expansion stage. Therefore in order to compare the theory with the observation, 
it is necessary to solve the evolutions of the cosmological perturbations on superhorizon 
scales using the concrete theoretical models such as the various inflation and the reheating 
scenarios. [3]] |1|] [3] [2]] P] [7]] |H]] P] Fortunately as for the evolutions of the cos- 
mological perturbations in the long wavelength limit, the exact solution is constructed in 
terms of the evolution of the corresponding locally homogeneous universe; as for the linear 
perturbations in the papers [lO]] [5]] [11]] [8]] and as for the full nonlinear perturbations 
in the papers [13] ] [M] ] P ] . The final form of the exact solutions of the evolutions of 
the cosmological perturbations in the long wavelength limit is established by the Kodama 
Hamazaki construction (KH construction), as for the linear perturbations in the papers [5] 
] [8]], and as for the full nonlinear perturbations in the paper P]. In the KH construc- 
tion, the physical quantities related with the exactly homogeneous universe, such as the 
scalar quantity perturbations, are given as the solutions of the evolution equations of the 
corresponding locally homogeneous universe and the physical quantities not related with 
the exactly homogeneous universe, such as the vector quantity perturbations, are given by 
solving the first order evolution equations, that is, the spatial components of the Einstein 
equations. It was shown that the second order evolution equations of the spatial unimod- 
ular metric including the information of the adiabatic decaying mode is exactly solvable. 
In the present paper, we use the KH construction. 

The general theory of relativity is a gauge theory. When we solve the equations of the 
general theory of relativity, the nondynamical gauge modes are contained in the solutions. 
Therefore in order to extract the dynamical modes only, it is desirable to write down the 
equations in terms of the gauge invariant variables only. In the linear perturbation theory, 
the program of the gauge invariant perturbation variables was first performed in the paper 
[T] ] , and was extended so that we can treat the multicomponent systems [TT] ] [TS] ] [IS] 
]. In the second order perturbation theory, the gauge invariant perturbation theory was 
constructed in the papers [36] ] [37] ] [15] ] . In our previous paper P ] , the full nonlinear 
perturbation theory in the leading order of the gradient expansion was constructed and 
several main definitions of the gauge invariant perturbation variables including the nonlin- 
ear Bardeen parameters [I] ] [IH] ] [12] ] [S] ] P ] were presented. In the present paper, in a 
more general way, definitions of all the types of the gauge invariant perturbation variables 
are constructed and it is manifestly shown that all the perturbation equations of the phys- 
ical laws such as the evolution equations, the constraint equations and the metric junction 
conditions can be written by using the gauge invariant perturbation variables which we de- 
fined only. By solving the equations of the gauge invariant formulation of the full nonlinear 
perturbations in the leading order of the gradient expansion which we formulated, we can 
extract the full nonlinear physically meaningful, dynamical information of the cosmological 
perturbations on superhorizon scales. 



In order to interpret the physics of the evolutionary behaviors of the cosmological 
perturbations, the Adiabatic/Entropic decomposition (A/E decomposition) of the cosmo- 
logical perturbations is efficient. The essence of the A/E decomposition is in defining 
the adiabatic perturbation variable and the entropic perturbation variables. Although 
the linear version of the A/E decomposition has been already established [T7| ] [H]], the 
satisfactory definitions of the adiabatic perturbation variable and the entropic perturba- 
tion variables in the nonlinear perturbation theory has not been completed yet. In the 
present paper, we give the definitions of the adiabatic perturbation variable and the en- 
tropic pertubation variables which can be used in the nonlinear perturbation theory by 
using the fact that the universal adiabatic growing mode always exists in the solutions in 
the nonlinear perturbation theory in the long wavelength limit. [5' ] That is, we call the 
perturbation variable which does not vanish for the universal adiabatic growing mode the 
adiabatic perturbation variable and we call the perturbation variable which vanishes for 
the universal adiabatic growing mode the entropic perturbation variable. In particular, the 
adiabatic/entropic perturbation variables which are defined under the p philosophy where 
the evolutions of the system are traced by choosing the energy density p as the evolution 
parameter, have desirable properties. All the perturbation variables in this set are con- 
tinuous across the metric junction hypersurface which is defined by p = const such as the 
slow rolling-oscillatory transitions of the scalar fields and the reheating transitions. The 
evolution equations of the perturbation variables in this set which can be derived by quite 
easy calculation have very simple expression. The adiabatic perturbation variable in this 
set is the well-known Bardeen parameter. P] [19]] [l2]] [M]] P] 

In the near future, more precise observations of CMB will be performed and it is 
expected that the information of the nonlinearity of the CMB fiuctuations will be obtained. 
Motivated by the observational advancement, the models which generate the significant 
nonlinearity characterized by the large non-Gaussianity parameter f^i [35| ] have been 
proposed; the inhomogeneous end of the infiation |28]] [29|], the modulated reheating [31] 
], the curvaton scenario j32]], the vacuum dominated infiation [53]] [SI]]- The former two 
cases are related with the metric junction hypersurface which cannot be defined by p = 
const and the large non-Gaussianities f^L are generated discontinuously on the transition 
hypersurface. In the latter two cases, the non-Gaussianity f^L grows continuously and 
becomes very large transiently. In the present paper, we present the exponent evaluation 
method which enables us to give the analytic order estimates of the non-Gaussianities 
/jvL in these models. We discuss that the mechanisms which generate the large non- 
Gaussianities /at^, in the latter two cases are common, although in the first case in the 
latter two cases the cosmological term does not exist while in the second case in the latter 
two cases the cosmological term exists. In the latter two cases, the entropic perturbation 
of the component which does not govern the cosmic energy density can trigger the growth 
of the Bardeen parameter Cn(p) and the non-Gaussianity Jmli when the scalar fields which 
support the entropic perturbation are very small, since the influences of these small scalar 
fields on the Bardeen parameter Cn(p) can become large. 

The rest of the present paper is organized as follows. In the section 2, we give the def- 
initions of all the types of the gauge invariant perturbation variables and show manifestly 
that in the long wavelength limit all the perturbation equations of all the physical laws 
derived by the general theory of relativity can be written in the gauge invariant manner. 



In the section 3, under the p philosophy, we complete the A/E decomposition of the full 
nonlinear perturbations by giving the definitions of the adiabatic perturbation variable 
and the entropic pertubation variables. In the section 4, as the application of the A/E 
decomposition based on the p philosophy formulated in the previous section, we investi- 
gate the evolutions of the cosmological perturbations in the universe where the growth of 
the adiabatic perturbation variable called the Bardeen parameter [T| ] [T9| ] [12] ] [H] ] |9] ] 
is induced by the entropic perturbation of the subdominant component. We evaluate the 
non-Gaussianity parameter f^i by the exponent evaluation method. We present the con- 
dition for which the non-Gaussianity fjsL becomes large in the models where fjsL changes 
continuously. 

§2 the manifestly gauge invariant formulation of the 
nonlinear cosmological perturbation theory in the 
leading order of the gradient expansion 

2.1 the evolution equations and the constraint equations 

We consider the Einstein equations G^^ = n'^T^,^ where k^ is expressed in terms of the 
Newtonian gravitational constant G as k^ = 8nG, using the 3 + 1 decomposition. [38|] [16] 
] P] The Greek indices p,^,- ■ ■ run from to 3 and the Latin indices i,j, ■ ■ ■ run from 1 
to 3. The metric tensor g^,y is expressed as 

^00 = -a' + /3fc/3^ (2.1) 

9o^ = A, (2.2) 

9ij = lij, (2.3) 

where a is the lapse and /3j is the shift vector. The index of /3j is raised by 7*-' which is the 
inverse matrix of jij. The spatial metric jij is factorized as 

lij = d^lij-. (2-4) 

where 7jj is the unimodular matrix whose inverse matrix is expressed as t'*-' and a is the 
scale factor. The energy momentum tensor of the total system T^^ is expressed as 

Ty^v = (p + P)uyUy + Pg^^v, (2.5) 

where p, P and m^ are the energy density, the pressure and the four velocity vector of the 
total system, respectively. Because of the normalization condition u^Uy, = —1, u^ can be 
parametrized as 

uo = -u'^ia^-P'iPk + Vk)}, (2.6) 

Uk = u\(3k + Vk), (2.7) 

where u^ = g^^u^ is given by 

u' = {a' - {f3k + Vk)if3' + v')y'/', (2.8) 



Vk is the three velocity of the total system and the index of v^ is raised by 7'-^. T^,^ is 
expressed as 

a 

where Ta^^, is the energy momentum tensor of the perfect fluid component a and Tsf_iu is 
the energy momentum tensor of all the scalar fields. Tq^j, is expressed by (12. 5p (12. 6 p (12. 7p 
(12. Sp where p, P, u^ and Vi are replaced with pa, Pq, Ua/i and faj, respectively. Ts^y is 
expressed by 

T5^. = Y. ^^^'i'adAa -\\Y. 9"'dp(j)ad.^a + 2U i g^,. (2.10) 

a \ a, ) 

As for the scalar fields, since we cannot decide to which component a each term of the 
potential U belongs, Ts^u cannot be decomposed into Ta^v In this way, the indices of the 
component A are divided into the perfect fluid indices a and the scalar field indices a. The 
energy momentum transfer vectors of the perfect fluid component a and the scalar field 
component a are expressed by 

Qa^ = Sad^(l)a, (2.12) 

where Qa and /^^ are the energy transfer and the momentum transfer of the perfect fluid 
component a, respectively and Sa is the source function of the scalar field component a. 
The energy momentum conservation gives 

$^Qa/. + $^Qa^ = 0. (2.13) 

a a 

As for the perfect fluid component a, V ^T^^ = Q^v gives the equation of motion of the 
perfect fluid component a. As for the scalar fields, V^Tg^ — ^^ Qav = can be expressed as 
the linear combination of d^cpa- By assuming that the each coefficient of di^cpa is separately 
vanishing, we can derive the phenomeno logical equation of motion of the scalar field (pa, 

00, - dU/d(l)a = Sa. 

Since we want to treat the cosmological perturbations on superhorizon scales, we put the 
gradient expansion assumptions by using the small parameter e characterizing the inverse 
of the long wavelength of the cosmological perturbations. Since the spatial scale of the 
inhomogeneity of all the physical quantities is of the order of 1/e, we assign di = 0(e). As 
for the metric, we assign goi = 0(e). For arbitrary vector fields V^ satisfying V^Vfj, = 0(1) 
including m^, u^^, we assume that Vi = 0(e). Therefore /3j, (3\ Vi, v\ Vai, f^ and fai 
are of the order of e. As for the velocity vector of the total system and the perfect fiuid 
component a, the leading order of the gradient expansion can be expressed by 

Mo = -a + 0(e2), (2.14) 

«^ = -{vi + (3i) + 0{e^), (2.15) 

a 

Uao = -a + 0{e^), (2.16) 

Uci = -iVa^ + P^)+Oie^). (2.17) 

a 



As for the momentum transfer vector of the perfect fluid component a, u^fafj. = gives 

/„0 = + O(e2). (2.18) 

We consider the gauge transformation laws of all the physical quantities. The gauge 
transformation laws are written in terms of the Lie derivative. The Lie derivatives of the 
quantity with upper index and the quantity with lower index are expressed by 

L{T)X'' = TPdpX^ -dpT^'XP, (2.19) 

L{T)X, = T%X, + d^TPXp. (2.20) 

The Lie derivative of the tensor field of an arbitrary rank is given by the above two defini- 
tions and the Leibniz rule. Because of the gradient expansion assumption, the infinitesimal 
coordinate transformation generating the Lie derivative T^dp satisfies T* = 0(e). Under 
the gradient expansion assumption, the Lie derivative of the scalar S is given by 

L{T)S = T^S + 0{e^), (2.21) 

and the Lie derivative of the vector Vp is given by 

L{T)Vo = rVo + TVo + 0(e2), (2.22) 

L{T)Vi = T^V, + diT°Vo + 0{e^). (2.23) 

Under the gradient expansion scheme, it is possible that the quantity which is not a scalar, 
for example 'jij, has the Lie derivative of the scalar type (I2.2ip . So we expand the definitions 
of the scalar field and the vector field as follows. 

Definition The physical quantity which has the Lie derivative (I2.2ip is called the scalar 
like object. The physical quantity which has the Lie derivative (I2.22p (I2.23P is called the 
vector like object. 

Following these definitions, the physical quantities such as a, jij, jij, 7*-', 7*-^ are the 
scalar like objects and the dp derivative of these quantities is the vector like object. We 
can demonstrate the following propositions easily. 

Proposition 1 For a scalar like object S, dpS is a vector like object. 
Please use L{T)dpA = dp{L{T)A} for an arbitrary quantity A. 

Proposition 2 For two arbitrary vector like objects Ap, Bp, 

#, A-^B. (2.24) 

-Do -Do 

are scalar like objects. 

Corollary For a scalar like object A, DfA where 

1 d 
A := -^ (2.25) 

act 

is also a scalar like object. 

Please notice that uq = —a + O(e^), and that dpA is a vector like object. 



Corollary For scalar like objects A, B, Di(A)B where 



A(A):=9,-^|, (2.26) 



is also a scalar like object. 

Please notice that df^A, d^B are vector like objects. 

Corollary For a scalar like object A, 



d,A+—{v, + Pi) (2.27) 



a" 



is also a scalar like object. 

Please notice that the above quantity can be written as diA — {A/uQ)ui where u^ is the 
velocity vector of the total system and that d^A is a vector like object. 

Proposition 3 In the background level, all the evolution equations and all the constraints 
can be expressed in the form that polynomials of the scalar like objects only are vanishing. 
As the proof, we write down the Einstein equations. As for the space-space components 
of the metric tensor, we use the matrix notation: M := (7ij), M~^ := (7*-'). H is the 
Hubble parameter defined by a/ a The Einstein equations G^^ = k^T^^ give the Hamiltonian 
constraint 

(-h) ='lp + l^tT(-MM~'-MM-'), (2.28) 

\a J 3 24 \a a J 

and the evolution equations 



(-) = -itr f-MM-'-MM-') - ^(p + P), (2.29) 



ia/iaAA^3£/l^.x_i^.^,^_,i^.^_^_ p^3^l 



a dt \a dt J a \a J a a 

and the momentum constraint 



3 ^ ^ J 

1 " 
'2 



M-'-dt{D,{a)M} 
a 



1 / 1 . 

+ -tr M-^DAa)M ■ M-^-M 



+2Dj{a) [ — ]- i^^h--Zj, (2.31) 



a J aH 

where Zi is the scalar like object defined by 



Zi:=dia + -^ivi + Pi). (2.32) 



n0„ - dU/d(j)a = Sa gives 

1 d /I 



ija \ H (ba dU 



adt \a dt J ' " a a ' dcpa 
As for the perfect fluid components, V^T^g — Qao and V^T^j = Qai give 

1 M 

Pa = -3 — (Pq + Pa) + <5a, 

a 



(2.33) 



and 



a 



■ 2, ^7 

'^ '^ TT '^^ 

ri 



(2.34) 



T^j'^a^i Ji 

an 



Di(a)Pa + hatt — Di(a) ( — 
a \a 



(2.35) 



where h^ := pa + Pa is the enthalpy of the fluid component a and Zai is the scalar like 
object defined by 

(2.36) 



a 



Zai ■= ditt + —{Vai + A)- 



a' 



Then we conclude that all the evolution equations and all the constraints can be written 
in terms of the scalar like objects only. The evolution equations of M (I2.30p can be solved 



as 



M = Ri exp 



io 



dt^R2 



(2.37) 



where Ri, R2 are the 3x3 time independent matrices depending on x: Ri is unimodu- 
lar symmetric, R2 is traceless and R1R2 is symmetric. [9]] By using( p37|) . the term in 
(I2:28|)(l2:29|) can be written as 



where 



-ti(-MM-^-MM-A =% 
4 \ a a a° 



cr ■= -MRl 



(2.38) 



(2.39) 



We consider the perturbation. We assume that the arbitrary background quantity A 
depends not only on (t, cc), but also on A which characterizes the perturbation. We can 
Taylor expand A around A = as 



•4(A=i)=5:^ 



1 dJ'A^X) 



fc=0 



k\ dX'' 



(2.40) 



A=0 



where A{X = 1) is a full nonlinear quantity. We can identify 

d^A(A) 



dX'' 



^(5M, 



(2.4i; 



A=0 



where S^A is the k-th order perturbation of A. The gauge transformation of the background 
quantity A is defined by 

A{X, /i) = exp [fxL{T)]A{X, /i = 0), (2.42) 

where L{T) is the Lie derivative generated by the infinitesimal displacement T := T^d^, 
A{X, fi = 0) is the quantity before the gauge transformation and A{X, /x = 1) is the quantity 
after the gauge transformation. This expression (I2.42p is a solution of the differential 
equation 

-^A = L{T)A, (2.43) 

djj, 

which we use instead of the solution 02.42p from now on. By differentiating 02.43P with 
respect to A, we get 

d dA fdT\ ^ , ,dA , , 

since not only the background quantity A but also the infinitesimal displacement generat- 
ing the Lie derivative T depends upon A. In general, the gauge transformation of the A 
derivative of A contain the Lie derivatives generated by d^T /dX^. But we can make a new 
quantity B by combining the A derivatives of the background quantity appropriately, so 
that 

^B = L{T)B (2.45) 

dfi 

which does not contain the Lie derivatives generated by d^T /dX^ (A; = 1, 2, ■ ■ ■ ) can hold. 
So we can put the definition as follows. 

Definition We call a quantity B which has the gauge transformation (I2.45P the back- 
ground like object. 

Any background like object is a gauge invariant quantity with respect to all the in- 
finitesimal gauge transformation satisfying T(A = 0) = 0. We can prove the following 
proposition. 

Proposition 4 Let A, B be the scalar like objects and the background like objects. Then 
D{A)B where 

D{A) -.= 4,- ^i#, (2.46) 

^ ^ dX dXAdV ^ ^ 

is also a scalar like object and a background like object. 

Corollary Under the assumptions in the previous proposition, D{A)"'B {n = 1,2,---) 
are also scalar like objects and background like objects. 
For the proofs, please see our previous paper. [9]] 



Proposition 5 All the perturbation equations of the evolution equations and the con- 
straints can be expressed in the form that the polynomials of the quantities which are the 
scalar like objects and the background like objects are vanishing. That is, all the per- 
turbation equations of the evolution equations and the constraints can be expressed in a 
manifestly gauge invariant manner. 

The perturbation equations can be obtained by operating D{S) where S is an arbitrary 
scalar like object on (Km (K29\) (lOOl) (lOTll I^M) (EMD (ESSD which are written in 
terms of the scalar like objects. Since a scalar like object operated D{S) on is a scalar like 
object and a background like object, the assertion of the proposition 5 follows. If we want 
to move the time derivative to the outermost position, you can use 



adt 



-^DiS) (^] i|, (2.47) 

a ^ ' \S adt ^ ' 



where [A, B\ := AB - BA. 



2.2 the junction conditions 

In the early universe, there exist periods when the equation of state changes quite rapidly, 
such as the slow rolling-oscillatory transition and the reheating transition [4 ]. As the 
zeroth order approximation, it is appropriate to treat these transitions by connecting two 
spacetimes which have different equations of state by the metric junction formalism [23]]. 
These transition hypersurfaces are defined by the particular equations; for the slow rolling- 
oscillatory transition, Hja = m, and for the reheating transition H/a = T where m, T are 
the mass, the decay constant of the scalar field, respectively. Motivated by the above point, 
we extend the metric junction theory across the spacelike hypersurface defined by C = 
where C is the scalar like object, within the framework of the full nonlinear perturbation 
theory in the leading order of the gradient expansion. In the appendix B, we formulate the 
metric junction in the linear perturbation theory in the long wavelength limit, and mention 
its consistency with the full nonlinear theory. 

We consider a 4 dimensional spacetime Ai and a 3 dimensional hypersurface S. S 
separates Ai into two region: Ai^ which is the future of S and A^_ which is the past of 
E. The hypersurface S is parametrized by the intrinsic coordinate y* (i = 1, 2, 3) as 

xl = t^+6Z±{y), (2.48) 

x'^ = y' + 6ZUy), (2.49) 

where x^ are spacetime coordinates in the region Ai±, respectively and tx is a constant 
common to A^j-. From now on, we omit index ±. The gauge transformation of 6Z, 6Z^ 
are given by 

L{T)6Z = -T°, (2.50) 

L{T)6Z' = -T\ (2.51) 

Proposition 6 Let A^, S^ be vector like objects. Then A^diS Z + {Aq / BQ)Bi is the scalar 
like object. 



Corollary 0j := adiSZ + {a/a)dia is a scalar like object. 

In the previous proposition, as A^, B^, please adopt m^, d^a, respectively. 

As the junction hypersurface, we adopt the hypersurface characterized by C = where 
C is a scalar like object. Then we get 



di6Z 



which yields 



a 



C 



diC 



Di{a)C. 



The normal vector n^ of the hypersurface S pointing from M._ to A^+ is given by 

-sgn(C') 



Un 



-gp-dpCd^C 



d,C, 



and the tangential vector ef on S are given by 






1,2,3). 



Then we get 



n^n'' 



-1, 



n^ef = 0. 



We define the intrinsic metric qij and the extrinsic curvature Kij of S by 

Kij := efe'^jVf.n^. 
In case of S defined by C = 0, we obtain 









As for the energy momentum tensor T^,^, we obtain 



T 

T 



«J 



n^^e'^Tp, = -{p + P)-{Z, - D,{C)a] + 0{e' 

(Jj 



We notice that q^j, K^j, Tnn, T^i and Tij can be written by the scalar like objects only, 
junction condition formulated by Israel [231] is given by 

[q,j\t = [K,j\t = [TrX- = [Tr.:\t = 0, 
where [Qjt := Q+ — Q-- In our notation, the above juction condition is written by 



[7 



»jj 



a 



[p\t = [{p + P){Z,-D.-.{C)a}]t = ^. 



2.52) 
2.53) 

2.54) 



2.55) 

2.56) 

2.57) 
2.58) 

2.59) 
2.60) 

2.61) 

2.62) 

2.63) 

The 

2.64) 
2.65) 



10 



Proposition 7 In the background level, the metric junction condition can be expressed 
in terms of the scalar like objects only. 

We consider the perturbation of the junction condition. As for the perturbation, the 
next proposition is essential. 

Proposition 8 Let the matching hypersurface be defined by C = where C is a scalar 
like object. For an arbitrary scalar like object S satisfying [S]t = 0, D{C)S, Di{C)S are 
continuous across the matching hypersurface: [D{C)S]t. = 0, [Di{C)S]t = 0. 

For the proof, please see the appendix A. By applying the above proposition finite 
times, we obtain the following corollary. 

Corollary Under the assumption presented by the previous proposition, [D{C)^S]t = 
for an arbitrary natural number n. 

As for M := (%•), M is solved as fl237l) . From fl235D . M in A^+ is given by 



M+ = Ri_ exp 



dt—R2+ / dt—R2 



(2.66) 



where i?i_ is -Ri in A^_, and i?2 in ■M.+ and R2 in Ai_ is the same R2+ = R2- ='■ -^2- 

As the junction, we consider the transition where the energy p^_ transfers into the 
energy pA+ which has the different equation of state from pa-- From fl2.65p . the energy 
momentum conservation 

[pA]t = 0, (2.67) 

[{pA + Pa){Zm - A(C)a}]l = 0, (2.68) 

must hold. In the above discussion, all the perturbation equations of the metric junction 
conditions are written in the form that the polynomials of the quantities which are the scalar 
like objects and the background like objects are vanishing, therefore all the perturbation 
equations of the metric junction conditions are gauge invariant. 

§3 choice of the independent gauge invariant variables 
based on the classification of the perturbation so- 
lutions into the adiabatic mode and the entropic 
modes 

3.1 the universal adiabatic growing mode 

All the evolution equations of the locally homogeneous universe are invariant under the 
transformation defined by 

a -^ aA, (3.1) 

R2 -^ i?2A^ (3.2) 

a — 7- a, (3.3) 

W ^ W, (3.4) 

11 



where a is the lapse function and W is an arbitrary scalar quantity. A is a time inde- 
pendent function A = K(x\ Taking the variation with respect to A considering that only 
A is dependent upon A gives the obvious perturbation solution. We call this solution the 
universal adiabatic growing mode. [5]] For an arbitrary scalar quantity S", the first order 
and the second order perturbation solutions of the universal adiabatic growing mode is 
written as 

D{a)S = -|^lnA, (3.5) 

-<")'^ --11^-41 (I) (1;")^ <-) 

where the first order expression is very familiar in large literature. We call the gauge 
invariant perturbation variable defined by 

US) ■.= D{Sr In a (3.7) 

the generalized Bardeen parameter induced by the scalar like object S. When we adopt an 
arbitrary scalar quantity W or H/a as the scalar like object S, the perturbation solution 
of the universal adiabatic growing mode is written as a time independent form: 

Jn 

US) = —lnA. (3.8) 

3.2 the adiabatic perturbation variable and the entropic pertur- 
bation variables 

In order to interpret the physics of the linear cosmological perturbations, the classification 
into the adiabatic perturbation and the entropic perturbations was often convenient. [17] 
] [18] ] [19] ] [2] ] Therefore the generalization of this classification into higher order pertur- 
bations are thought to be useful. So we define the adiabatic perturbation variable and the 
entropic perturbation variables in the higher order perturbation theory. 

Definition We call the perturbation variable which does not vanish for the universal 
adiabatic growing mode the adiabatic perturbation variable. We call the perturbation 
variable which vanishes for the universal adiabatic growing mode the entropic perturbation 
variable. 

We will present the examples of the adiabatic and the entropic perturbation variables. 
We assume that S, Si {i = 1, 2) are the scalar like objects such as W, W/a, H/a where 
W is an arbitrary scalar variable. The generalized Bardeen parameter (n{S) and D{a)"'S 
are adiabatic perturbation variables and Cn('S'i) — Cn(5'2) and D(5'i)"'5'2 are entropic per- 
turbation variables. 

3.3 the N philosophy and the p philosophy 

We call the expressions representing the physical quantities at the final time in terms of 
those at the initial time the S formulas. [8]] Our final purpose is to construct the S formulas 
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of the adiabatic perturbation variable such as the Bardeen parameter Cn(p) •= D(p)"'\na. 
In the previous subsection, it was shown that the Bardeen parameter Cn(p) is time indepen- 
dent for the universal adiabatic growing mode. Therefore we expect that the formulation 
in which the difference between the Bardeen parameter at the final time and that at the 
initial time can be expressed in terms of the entropic perturbation variables may exist. In 
this subsection, we choose the appropriate set of the entropic perturbation variables and we 
construct the formulation in which the time change of the Bardeen parameter is brought 
about by the evolutions of the set of these entropic perturbation variables. 

Until now in order to understand the evolutions of linear cosmological perturbations in 
the universe governed by the multiple component energy densities, the decomposition of 
the perturbations into the adiabatic component and the entropic components has already 
been performed. [IT] ] [H] ] [T9| ] [2] ] In the nonlinear perturbation theory, the following set 
of perturbation variables was adopted: as the adiabatic perturbation variable, the Bardeen 
parameter Cn{p), and as the entropic perturbation variables, the difference between the 
generalized Bardeen parameters induced by the energy densities of the different compo- 
nents Sn{pA, Pb) '■= CnipA) — Cnips) whcre the subscripts A, B represent the different 
components. [20j ] Since all the perturbation variables in this formulation [20] ] are based 
on perturbations of the logarithm of the scale factor A^ := In a, we call this formulation 
the A^ philosophy. But in the A^ philosophy, it is difficult to write down the evolution 
equations in terms of the set of variables (nip), Sn{pA, Pb) in the closed form. From now 
on, we often consider the matching of the metric across the matching hypersurface defined 
by p = const. Across such matching hypersurface, the perturbation variables in the A^ 
philosophy, Sn{pA, Pb) jump by finite values. 

In order to solve the defects in the A^ philosophy, we propose the new set of the per- 
turbation variables. We choose 



-j l„a. d[^-) s., (3.9) 

where sa '■= Pa/p, as the adiabatic perturbation variable and the entropic perturbation 
variables, respectively, sa satisfies J2a^a ~ ^- -^^ ^°^ ^^^ ^^^ ^^^ °f ^^^ perturbation 
variables, no finite jumps do not exist across the slow rolling-oscillatory transition H/a = m 
where m is the mass of the scalar field and across the reheating transition H/a = T where 
r is the decay constant of the scalar field. In order to avoid the calculational complexity, 
we assume that c/j = 0, since we can neglect the second term of the right hand side of the 
Hamiltonian constraint fl2.28p with f l2.38p because of the rapid growth of the scale factor 
a during the inflationary expansion of the universe. In the condition c/j = 0, the matching 
conditions of the slow rolling-oscillatory transition, of the reheating transition are reduced 
into p — "iw? / H? = 0, p — SF^/k^ = 0, respectively. Under the simplification of cr = 0, the 
set of the perturbation variables which we adopted in (13. 9p is reduced into 

D(p)"lna, /^(p)"sA. (3.10) 

Since these perturbation variables are continuous across the matching hypersurface defined 
by p = const, we only have to concentrate on solving the evolution equations of these 
perturbation variables. Since all these variables are defined by D{p), we call the use of 
these variables presented in f l3.10p the p philosophy. 
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We will give the evolution equations of the perturbation variables fl3.10p . For simplicity, 
we assume that the multiple components do not interact and that pa obeys dp^/dN = 
—QaPa where qa will be called the g factor from now on. When the a component is 
the perfect fluid with Wa '■= Pa/ Pa, its g factor is given by ga = 3(1 + Wa)- When the 
a component is the slow rolling massive scalar field with mass rria, its g factor is given 
by ga = 2ml/K,'^p. Since it was shown that the oscillatory massive scalar field can be 
approximated by the perfect fluid with w^ = [3] ] P ] P ] P ] , we can use the perfect 
fluid with ga = 3 instead of the oscillatory massive scalar field. The evolution equations of 
N ■.= \na, sa are given by 

^N = -1, (3.11) 

dp ps 

d 1 / gASA\ ,„ ION 

—Sa = - -saH , (3.12) 

dp p \ s / 

where s := YIbSb^b- We choose the total energy density p as the evolution parameter 
instead of the cosmic time t, and the right hand sides of (13. lip . (I3.12p are written by p, sa 
only. The evolution equations of the perturbation variables in the p philosophy are given 
by operating D{p) finite times on (13. lip . (I3.12p . In this case, it is important to notice that 
D{p) and d/dp are commutative since d/dt and d/d\ are commutative. 

In the p philosophy, all the perturbation variables are continuous across the matching 
hypersurface defined by p = const because of the proposition 8, and we can easily derive 
the evolution equations of the perturbation variables. Since the p philosophy is superior to 
the N philosophy because of the above two reasons, we will adopt the p philosophy from 
now on. 

§4 the non-Gaussianities of the nonlinear cosmologi- 
cal perturbations 

In this section, we discuss the non-Gaussianities generated in several cosmological models. 
The non-Gaussianities are measured by the f^L parameter. [35] ] It is assumed that the 
logarithm of the scale factor A^ := Ina is given by the function of the energy density p 
as the evolution parameter and of the solution constants. We only consider the models 
where the origins of the cosmological perturbations are in the quantum fluctuations of the 
scalar fields (pa in the inflationary universe. The statistical mean values of the perturbation 
amplitudes are given by 

'dct>a{^)dM^\\ , 

where (/>a(0) is the expectation value of the scalar field (pa at the first horizon crossing and 
H is the Hubble parameter at the first horizon crossing. In this case, the solution constants 
are given by the set of the expectation values of the scalar fields at the first horizon crossing 
{0a (0)}. In this case, using the logarithm of the scale factor A^ = N{p, 0i(O), 02(0), ■■■ ) 
the non-Gaussianity parameter /tvl is defined by 
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d(f)a{0) d(j)a{0)d(f)b{0) 

Na, Nab are defined as the coefficients given when we expand the gauge invariant adi- 
abatic perturbation variables D{p) In a, D(p)^lna with respect to d(f)a{0)/d^, respectively; 

fl(rtl„a = 5:iV„f^, (4.4) 

a 

D[p) lna = 2^iVa6— 77 ^7— • (4.5) 



ab 



We assumed that the more than second order perturbations of (pa at the initial time are 
all vanishing; (i"0a(O)/rfA" = (n > 2). 

Since the cosmological perturbations have the origin in the scalar field fluctuations 
in the de Sitter stage, the deviations of the spectral indices of the Bardeen parameter 
Ci(p) from the scale invariance (iln<< Ci(p) »/d\Y].k are suppressed by the slow rolling 
parameter. Since the first horizon crossing is defined by the relation 

k = aH = e^' ^p{QY'\ (4.6) 

where A^*, p(0) is the logarithm of the scale factor, energy density at the first horizon 
crossing time, respectively, we obtain 

The slow rolling phase is characterized by the smallness of the g factors of the scalar 
fields 0a(O) whose size is bounded by 5$ a small constant characterizing the slow rolling 
of the scalar fields; \ga\ < Ss where ga is defined by the evolution equations of the energy 
densities of the scalar field (pa', Pa at the first horizon crossing: dpa{0)/dN^ = —gaPaiS^)- In 
all the cases which we consider, the following evaluations hold: (9 In << Ci (p) »/dpa{0) ~ 
l/pa(0). By using the above properties, we can conclude that the Bardeen parameter in 
the first order perturbation theory Ci(p) has almost complete scale invariance: 

ln«Ci(p) »~(5s. (4.8) 



d\nk 



Then we concentrate on the non-Gaussianity parameter f^L from now on. 

Except for the cases where the large f^i is generated discontinuously on the transition 
hj^ersurface such as the inhomogeneous end of the infiation [28] ] [29] ] and the modulated 
reheating [30]] [31]], different two cases have been discussed. One is the curvaton scenario 
[32] ] and the other is the vacuum dominated two scalar fields [33] ] [M] ] • We discuss 
that the mechanism which generates the large f^i continuously in the above two different 
cases can be explained from the three common viewpoints which will be presented in the 
subsection 2 of this section. In this section, we use two strong methods; the exponent 
evaluation method and the p philosophy. By the exponent evaluation method, it becomes 
possible to evaluate the order of /at^ analytically in the wide ranges of parameters. In the p 
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philosophy, the time evolutions of the above two systems are traced using the logarithm of 
the scale factor N as the adiabatic independent variable, S2 which implies the ratio of the 
energy density of the component which does not govern the energy density of the universe 
as the entropic independent variable, and the energy density p as the evolution parameter. 
The p philosophy makes the instant when j'nl grows large and the length of the period 
when the large /nl continues clear. 



4.1 the exponent evaluation method 

In many papers, the calculations of the non-Gaussianity parameter Jnl were performed: for 
the inhomogeneous end of the inflation [28]] [291], for the modulated reheating [31]], for the 
curvaton [32]] and for the vacuum dominated inflation represented by the hybrid inflation 
[33] ] [M] ] . In many papers so far, the calculations of f^L were performed numerically and 
in the extreme situations where only one factor is concerned the analytic formulas of Jnl 
were given. In this subsection, we present the exponent evaluation method by which it 
becomes possible to give the analytic order estimates of fiyL in the wide range including 
cases where more than two factors are concerned. 

We explain the exponent evaluation method by adopting the inhomogeneous end of the 
inflation [28,] [29]] as an example. We consider the two scalar fields 0i, 02 governed by the 
vacuum dominated potential given by 



U 



2 I 



>l 



(4.9) 



a=l 



where rja {a = 1,2) are negative and Uq is the large constant compared with the terms 
quadratic with the scalar fields. Under the approximation where the vacuum energy Uq is 
the dominant contribution of the energy density p and where the scalar fields 0^ (a = 1, 2) 
are slow rolling on the potential, the evolutions of 0^ (a = 1, 2) are given by 



6„(0)exp 



Va 



K^a 



■N 



(4.10) 



For simplicity, rja (a = 1, 2) are assumed to be a independent: rja = rj. We assume that the 
inflation ends in the bifurcation set defined by 



J2^a<Pl 



a 



(4.11) 



On the bifurcation set, the waterfall field which interacts with the inflatons 0^ (a = 1,2) 
gets the negative mass and the large vacuum energy ~ Uq is transferred into the oscillation 
energy of the waterfall field and into the radiation energy which interacts with the waterfall 
field. Also in the situation where along the curve defined by (14. lip the deep ditch of the 
potential exists, the inflation ends on the curve (14. lip . The logarithm of the scale factor 
A*" at which the scalar fields reach the bifurcation set (14. lip is given by 



N 



kHIq 
2rj 



In 



^E^«<^'(o) 



(4.12) 
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In order that the inflation can solve the horizon problem, we assume that a := f^V^jIri 
10^. The non-Gaussianity parameter /jvl is calculated as 

fNL = -\^-A=---{v-^). (4.13) 



a \ 2Ao I a 



where 

2 

^n:=Y.^im\ (4.14) 

From now on, we neglect numerical factors of order unity without mentioning. By setting 

7r- := — = 10 , 0r := — 7T7T = 10 , 4.15 

71 0i(O) 

we obtain 

The exponent evaluation method gives 

k<\l P=l, (4.17) 

^/ < A; < / j9 = lO^'^^^', (4.18) 

l<k<2l p=lO-^^'^^, (4.19) 

2/ < A; p = 1. (4.20) 

As an example of application of the exponent evaluation method, we adopt I < k < 21. 
Since I < k < 21, we obtain 

1 + 10^-2' ~ 1, (4.21) 

1 + lO^'^-^i _ io3fc-2i^ ^4 22) 

1 + 102^^-2' ~ lO'^^-^i^ (423) 



and get 



1 1 C\3k—2l 



p takes the maximum at fc = /: p = lO', then /tv^, = 10'~^. In the exponent evalua- 
tion method, we take only the term which has the largest exponent in the polynomial 
constructed by several 10*^ type terms. 

As the second example of application of the exponent evaluation method, we consider 
the modulated reheating. [30]] [31]] The scalar field 0i on the potential U = m'^(j)l/2 causes 
the chaotic infiation: 

N = ^{<Pl{0)-<Pl}. (4.25) 

When H^ = m? where H is the Hubble parameter, the scalar field 0i begins to oscillate 
and behaves like a dust fiuid [3] ] [6] ] [7] ] [8] ] : 

^=^^?(0)-^-3>n(^). (4.26) 
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When H^ = F^ where F is the decay constant of the scalar field 0i, the scalar field oscillation 
is transformed into the radiation fiuid: 

iV^^^f(0)-f-il„^-i,n(|;g). (4.27, 

In the modulated reheating, we consider that the decay constant of the first scalar field 0i; 
F is the function of the second scalar field 02 as F = 0^02(0) where a^ is a constant and 
n is an integer. In this case, A^ takes the form as 

iV = fi:202(o) + ln02(O), (4.28) 

up to the p dependent part which does not contribute the Bardeen parameters Cn(p)- In 
order that the 0i infiation can solve the horizon problem, we assume that 0i(O) = IO/k. 
We assume that the second scalar field 02(0) takes a small value as 02(0) = 10^''/k. The 
non-Gaussianity parameter /tvl is calculated as 

10^ - 10^' 
f^^ = (102 + 10202- (4-29) 

The exponent evaluation method gives 

1 , 

1 <l 

When the second scalar field 02 takes a very small value, a significant non-Gaussianity J'nl 
is generated. 

4.2 the non-Gaussianities induced by the entropic perturbation 
of the component which does not govern the cosmic energy 
density 

In this subsection, we investigate the mechanism which triggers the large Jnl in the dif- 
ferent two models 1, 2 where the large Jnl can be generated continuously. The model 1 
is the radiation-dust system and the non-Gaussianity /jvl in the model 1 was investigated 
in the paper [32| ] in the context of the curvaton scenario. The model 2 is the vacuum 
dominated two scalar fields and the non-Gaussianity f^i in the model 2 was investigated 
in the papers [33] ] [311 ] in the context of the hybrid infiation. Although the model 1 and 
the model 2 are quite different apparently, we consider the mechanisms which generate the 
large /atl's in the model 1 and in the model 2 are completely the same. It is assumed that 
the infiation sufficient to solve the horizon problem A^ ~ 10^ is brought about by the first 
scalar field 0i. Under this assumption, the common points are summarized in the following 
three points: 

A. The expectation value of the second scalar field 02; 02(0) is very small compared with 
that of the first scalar field 0i; 0i(O) at the first horizon crossing. 
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inl = lo~^ 


(4.30) 


fNL = -lO-'^'-'\ 


(4.31) 


Inl = -1, 


(4.32) 



B. The component which originates from the second scalar field 02(0) does not govern 
the cosmic energy density p. 

C. The g factor of the component originating from the second scalar field 02(0); g2 is 
smaller than the g factor of the component governing the cosmic energy density p; 

The condition A guarantees that the contribution from the second scalar field 02(0) to 
the Bardeen parameters Cn(p) (^ = 1;2,---) is large. For example, we assume that A^ 
can be written as the sum of the 0i(O) dependent parts; A^ = Ylii fi{4>i{^)) ^^^ that each 
fi is the power or the logarithm of 0i(O). Then the contribution from 02(0) to Ci(p) 
is proportional to dN/d(j)2{0) ~ /2(02(O))/02(O) which is quite large when 02(0) is very 
small even if the contribution from 02(0) to A^; /2(02(O)) is quite small. Owing to the 
condition A, it becomes possible that the entropic perturbation of the component which 
does not govern the cosmic energy density p brings about significant contributions to the 
Bardeen parameters Cn(p) {n = 1,2, ■■■) without contributing A^. From the structure 
of the evolution equation of A^ := Ina fl3.1ip and that of sa fl3.12p . the growth of the 
Bardeen parameters Cn(p) '■= D{p)'"'N are governed by the entropic perturbation D{p)^SA 
of the component A whose ratio of the energy density sa is very small |sa| <^ 1. In 
addition, when \sa\ <^ 1, the entropy perturbation D[p)'^SA can grow or decrease rather 
rapidly. Therefore the condition B requires the existence of such component. Owing to the 
condition C; ^2 < gi, the contribution of S2 to A^ is monotonically increasing for the time 
evolution when the evolution parameter p becomes decreasing, therefore the contribution 
from D{pYs2 to the Bardeen parameter Cn(p) := D{p)"'N is also increasing. 

Since g2 < gi from the condition C, the ratio of the energy density of the second 
component S2 grows compared with that of the first component Si. Since the second 
component S2 begins to dominate the cosmic energy density p soon, the period when the 
condition B is satisfied is only the early period of time from the beginning. So the large 
f]\[L is realized only transiently in this period. 

Both the model 1 and the model 2 share the above three properties. The model 1 
is treated as the case without the cosmological term in section 4.2.1 and the model 2 is 
treated as the case with the cosmological term in section 4.2.2. 

In the rest of this subsection 4.2, we analyse the concrete physical systems using 
D{p)^\na, D{p)"-sa as the independent perturbation variables and p as the evolution pa- 
rameter. This scheme which we proposed as the p philosophy in subsection 3.3 is supported 
by the results of sections 2, 3. Subsection 2.1 guarantees that these independent perturba- 
tion variables D(p)"lna, D{p)"'SA are gauge invariant perturbation variables. Subsection 
3.2 states that D(p)"lna, D{p)"'SA can be regarded as the adiabatic perturbation, the 
entropic perturbations in the higher order perturbation, respectively. This A/E interpre- 
tation is useful when we interpret the time evolutions of the concrete physical systems. 
Since D[p)p = 0, that is -D(p) can be interpreted as the partial derivative with respect to 
A with p fixed, the p dependences of In a, sa are directly refiected to the p dependences 
of D{p)"'lna, D{p)'"'Sa- This fact makes the calculations and the interpretations of the 
time evolutions of the perturbation variables transparent. By the p philosophy which we 
explained in the above, in the rest of this subsection 4.2 we clarify that the entropy per- 
turbation D(p)"s2 supported by the energetically subdominant component p2 makes the 
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adiabatic perturbation D(p)"-\na and the non-Gaussianity /at^ grow considerably under 
the conditions that the g factor of this energetically subdominant component p2 is smaller 
than the g factor of the energetically dominant component pi and that the subdominant 
component p2 is supported by the extremely small scalar field expectation value. 

4.2.1 the case without the cosmological term 

We consider the two component system. We assume that the g factor of the component 
Pa (sa) is gA- Assuming that |s2| ^ 1 and linearizing f l3.lipfl3.12p with respect to S2, we 
obtain 

^-.2 = -(-l + '-^]s,, (4.33) 



dp p \ gi 

^N = -— + — f-l + ^~).2, (4.34) 

dp P9i P9i V 91. 



whose solution is given by 



., . 4^(^V'^"". (4.35) 



Pi(l) Vpi(i 



N = A^(l)--ln^ + -4^(^TT) , (4.36) 



-1+S2/31 



91 Pi(l) ^iPi(l) Vpi(1) 

where the g factors are assumed to be constant and X{1) implies the physical quantity X 
at an initial time. When g2 < gi, S2 and the contribution to A^ from the P2(l) dependent 
term increase and they are not bounded for the time evolution p — )■ 0. When g2 > gi, S2 
and the contribution to A^ from the P2(l) dependent term decrease for the time evolution. 
When g2 < gi, the ratio of the energy density S2 increases and reaches almost unity. In 
this case, the evolution equations (13. lip f l3.12p give 

".s, = l(-l + ii)s,, (4.37) 



dp p \ g2 

^N = -J- + J_(-l + «l).,. (4.38) 

dp P92 P92 V 92 J 

by linearizing (13. lip fl3.12p with respect to si assuming that |si| ^ 1. The solution is given 
by 

, . EM(^]-'^''"\ (4.39, 



P2(l) VP2(1) 

^..^W_l,„^,i^(^)-'^"'". (4.40, 

92 P2{1) 92p2{l) \P2{1)J 

When g2 < gi, si and the contribution to A^ from the pi(l) dependent term decrease for 
the time evolution p — )■ 0. 

In the above, the independent variables In a, sa are written as functions of p as the evo- 
lution parameter. By operating D{p) derivatives on expressions of In a, Sa, we can obtain 
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the expressions of the A/E perturbation variables D(p)"lna, D{p)'^SA in the form of the 
functions of p. Please notice that the D{p) derivative can be regarded as the partial deriva- 
tive with respect to A with the evolution parameter p fixed. Therefore the p dependences 
of In a, sa completely corresponds with the p dependences of D(p)"lna, D{p)^SA- In the 
above calculations and discussions, we can see that the entropic perturbation D{p)"'S2 of 
the energetically subdominant component p2 with g factor smaller than the g factor of the 
energetically dominant component pi makes the adiabatic perturbation variable D(p)"lna 
grow transiently while the energy density ratio S2 is increasing. 

When we calculate the non-Gaussianity /jvl, the expressions of A^'^, Nab are necessary. 
Na, Nab are defined as the coefficients of the gauge invariant adiabatic perturbation vari- 
ables D{p)\iaa, D{pY In a given when D{p)\na, D{p)'^\na are expanded with respect to 
d(l)a{0)/dX, respectively. The subscript a implies the d/d(f)a{0) derivative with p fixed. The 
p philosophy proposed in subsection 3.3 makes the calculations and the interpretations of 
the time evolutions of the non-Gaussianity J'nl as well as the A/E perturbation variables 
D{p)"'lna, D{p)"'sa more simple and more transparent. 

We apply the above results to the following concrete situation. Before A^ = A^(l), pi 
causes the chaotic inflation A^(l) ~ 10^, and at A^ = A^(l) decays into the radiation. After 
A^ = A^(l), pi is radiation [gi = 4). After A^ = A^(l), P2 is still in the slow rolling phase. 
Because we assume that 1712 ^ mi, 02 hardly moves from the initial value 02(0) {g2 = 0). 
At k'^p/S = m\, the slow rolling phase of p2 ends and begins to oscillate. After K^p/3 = m\, 
P2 behaves like dust fluid ((72 = 3). [3] ] P ] [7] ] [8] ] Then applying fl4.36p to the above 
situation yields 

Ar = A^ 1 +-h^-i^ + -mi0^ ^ — , o, , , ,,;, . 4.41 

^^4 p 8 ^^^^'\'iml) pi/4 (1 _ k202(O)/6)3/4 ^ ' 

This solution is simplified into the model given by 

Ar = fi;202(o) + fi;2^(p)02(o), (4.42) 

up to the p dependent part which is not related with the Bardeen parameters Cn{p), where 
a{p) is a function of p and increases for the time evolution p — )■ 0. All the numerical 
coefficients of order unity are dropped without mentioning from now on. We assume that 
(pi causes the inflation enough to solve the horizon problem and that 02 is a small field 
enough to contribute to the Bardeen parameters Cn(p) sufficiently: 

0i(O) = -10, 02(0) = -10-', (4.43) 

K K 

where / is a positive number. a(p) is written by a{p) = 10'^ where k increases for the 
time evolution p — t- 0. Since we adopt the approximation that p2 (^2) does not govern the 
cosmic energy density p, we obtain k < 21. The non-Gaussianity parameter /^r^ of the 
model fl4.42p is calculated as 

1 1 + l0^k-2l-2 

fNL = -r, ,, , ,^,,.,,_,,, - (4.44) 



102 (1 + 102fe-2«- 



21 



The exponent evaluation method yields 

2(/ + l) 



k< 

2{l + l) 



<k <l + l 



Inl = 10-^ 

1 Q3k-2l-4 



L 



NL 



l + l <k <2l 



L 



NL 
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-k+2l 



(4.45) 

(4.46) 
(4.47) 



For 2(l + l)/3 < k < l + l, f^i increases and reaches the maximum f^i = 10'~^ at k = l + l. 
For I + 1 < k < 21, fj^L decreases and reaches Jnl = 1 at k = 21. The non-Gaussianity 
parameter f^i takes a large value transiently. So if we want to obtain the large /tvl from 
the present observation, we need p2 to decay into radiation at k = I + 1. Next we consider 
the period when the second component p2 gets to dominate the cosmic energy density p. 
In this period, the first component pi is subdominant. So we can use (14.401) and get 



N = N{l) + -\n 



1 
+ 3 Vrni0i(O) 



m|0|(O) f«:Vi(l)f/' 
2p 

4/3 



3ml 



1 



;i - /t20i(o)/6)3/4 



3ml 



K^<Pl{0) 



6 



.1/3 



(4.48) 



which is simplified into the model (I4.28P for small p. Then we can obtain the evaluation 
of /jvL f O0|l ( lOT])fH:32|) . If we want /nl of order of unity, we need / > 1. 

There is a case where the contribution to A^ from the second component P2(l) dependent 
term increases but is bounded in spite of g2 < Qi. In this case, the non-Gaussianity 
parameter /tvl cannot grow into a significant value. We consider the case where (pa [a- = 
1, 2) with mass rria {m\ < mf) are in the slow rolling phase. In this case, Qa = 2m^/fi;^p. 
Unlike the previous case, ga depends on the cosmic energy density p. The ratio of the 
energy density of the second component S2 and the logarithm of the scale factor N are 
given by 



S2 



N 



P2(0) / p 



p(o) Vp(o) 



— l+mj/mj 



K 



2m{ 



(p-p(o)) + 



-mf + mi 



2m{ 



mr. 



P2(0) 



P(0) 



m\lm\ 



(4.49) 
(4.50) 



where in p(0), Pa(0) implies the first horizon crossing time and p(0) = pi(0) + p2(0) 
where Pa(0) = ?7i^0^(O)/2. From the above expression of A^, we can verify that the p2(0) 
dependent term is bounded and suppressed by i^^\{^\ therefore /atl is suppressed by 
10-2. 

4.2.2 the case with the cosmological term 

In cases where the cosmological term Vq exists, we use a := p — Uq as the evolution 
parameter. In these cases, the evolution equations corresponding to (13.111) . (I3.12p are the 
same evolution equations (13. lip . (13.121) but all p's are replaced with cr's. sa is defined by 
sa '■= Pa/o' and satisfies J2a^a = 1- When |s2| <^ 1, linearizing with respect to S2 gives 
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fl4.33p . f l4.34p but all p's are replaced with a's. In the same way, when |si| ^ 1, linearizing 
with respect to Si gives fl4.37p . (14.381) but all p's are replaced with cr's. 

We consider two scalar fields 0^ (a = 1,2) which move on the potential given by 



P = f/o + ^ Pa, Pa= ^Va(l)l- (4.51) 

a=l 

In this case, the g factor of the scalar field (j)a (a = 1, 2) is given by ga = 2ria/i^'^p = 
'^Va/i^^Uo. In order that the infiation can solve the horizon problem, we assume that 
ga ~ 10~^. We assume that rji > 772. First the energy of the scalar field 0i, pi dominate 
a and next the energy of the scalar field 02, P2 grows and gets to dominate a. In the first 
period where pi dominates a, N is given by 

2r,i pi(0) 2r,i pdO) \ <r J 



where 



Pa(0) := lva<Pl{0). (4.53) 



In the second period where p2 dominates a, N is given by 

i,^^^JEn,'.:f^EM(EMV-'"\ (4.54) 

2rh P2(0) 2ri2 p.JO) \ <r J ^ ' 

which can be obtained when all the subscripts 1, 2 are exchanged in the expression of A^ 
in the first period (I4.52p . 

We consider the case that 771 > 0, 772 < 0. In this case, in the first period a is positive, 
decreases for the time evolution, and then the p2(0) dependent term in A^ fl4.52p grows 
boundlessly. In the first period, A^ can be simplified into the model given by 

N = 10'^ln(j)i{0)-K^a{a)^l{0), (4.55) 

up to the a dependent part which is not related with the Bardeen parameters (n{p), which 
yields the non-Gaussianity parameter fisfL given by 

fNL = (-10'^ - «^V(c^)0^(O)) / (10'^ + ^'c^\^)4{0?j . (4.56) 

We put 

0i(O) = -lO'", 02(0) = -10-', a{a) = 10^ (4.57) 

K K 

where k grows for the time evolution and k < 2/+2 for the condition that p2 is subdominant. 
The non-Gaussianity parameter f^L is written by 

1 rj6— 4m i r\3k—2l 

J^^ ~ QQ4-2m _|_ lQ2k-2l\2' (^4.08J 
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The exponent evaluation method gives 



k < ki 












Inl = 


-lo-^ 


(4.59) 


ki < k < k2 












Inl = 


_iQ3fc-2i-8+4m 


(4.60) 


k2<k <k:i 












Inl = 


_lQ-k+2l^ 


(4.61) 


fci := 2 - 


4 


+ 


1-^ 


k,: 


:=2- 


- m + 1, 


fcs := 2/ + 2. 


(4.62) 



where 



For ki < k < k2, the absolute value of /jv^ grows and reaches the maximum /jvl = 10"^+"*+' 
at fc = k2. For k2 < k < k^, the absolute value of f^L decreases and reaches /tvl = 10^^ 
at A; = k^. In the second period, a is negative, —a grows, and the pi(0) dependent term in 
A^ (I4.54P decays. A^ in the second period (I4.54p can be simplified into the model given by 

N = -10^ In 02(0) + K^/3{a)(j)l{0), (4.63) 

up to the a dependent part which is not related with the Bardeen parameters Cn(p)- Putting 

0i(O) = ilO'", 02(0) = -10-\ /3{a) = 10-^ (4.64) 

K K 

where k increases for the time evolution and satisfies k > 2m — 2 for the condition that pi 
is subdominant, we obtain 

1 rj6+4i _j_ 1 ri—3k+2m 



{IQ21+4: _j_ lQ-2k+2m\2 
^ 10"^ (4.65) 

We consider the case r/i, r/2 < {rji > 772). In the first period, a is negative, —a increases, 
and the P2(0) dependent term in A^ (I4.52p grows boundlessly. In the first period, A^ (I4.52p 
can be simplified into the model given by 

N = -10^ In 01 (0) - k'^ a{a) (pliO), (4.66) 

up to the a dependent part which is not related with the Bardeen parameters Cn{p), which 
yields the non-Gaussianity parameter fj^fL given by 

Inl = (lO6^-«:V(a)0^(O))/(^lO^^ + /.V(a)0^(O)' 



^Q6+4m _ ^Q 



3k-2l 



ng2m+4 _|_ ]^g2fc-2n2' 



(4.67) 



by putting 



^(0) = -10^"^, 02(0) = -10-', a{a) = 10^ (4.68) 



K K 
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where in the first period pi is dominant m < I and k increases for the time evolution, 
k < 21 + 2 from the condition that p2 is subdominant. The exponent evaluation method 
gives 



k < ki 










Jnl = 10-^ 


(4.69) 


ki < k < k2 










/^^ = _io3fc-2'-4— 8^ 


(4.70) 


k2<k <k3 










fNL = -10-^+2', 


(4.71) 


2 
h := -^l + 


4 

-m 

3 


+ 2, 


k2: 


:=l + 


m + 2, A;3:=2/ + 2. 


(4.72) 



where 



For ki < k < k2, the absolute value of f^L grows, reaches the maximum f^L = — lO'"™'"^ 
at A; = k2. For k2 < k < /cs, the absolute value of /nl decreases, reaches /nl = — 10~^ at 
k = ks. In the second period, a is negative, —a grows and the pi(0) dependent term in A^ 
(I4.54p decays. N in the second period (I4.54p can be simplified into the model given by 

N = -10^ In 02(0) - ^2/3(^)0^(0), (4.73) 

up to the a dependent part which is not related with the Bardeen parameters Cn(p), which 
yields 



106+4« _ 10- 



3k— 2m 



(IQ21+4: _|_ 1Q-2k-2m\2 

10-2, (4.74) 



putting 



0i(O) = ilO-'", 02(0) = -10-', /3(a) = 10-^ (4.75) 

K K 

where k grows for the time evolution, and k > —2m — 2 from the condition that pi is 
subdominant. The above results about two scalar fields in the vacuum domination can be 
investigated by the method of the r function. P] The same results given in this subsection 
are reproduced using the r function in the Appendix C. 

The mechanism which produces the large /]vl depends on the fact that gi > g2 and 
that Si is dominant in the first time and that S2 begins to govern the cosmic energy density 
p gradually. So when si is the radiation gi = 4, and S2 is the scalar field with the negative 
mass g2 = 2ri/K,'^Uo {t] < 0), the non-Gaussianity Jnl can grow transiently, because the 
third term which depends on p2(0) 

can grow boundlessly. 
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§5 Discussion 

In this paper, the first half, that is sections 2, 3 is devoted to the general considerations 
about the gauge invariant nonlinear cosmological perturbation theory on superhorizon 
scales and the latter half, that is section 4 is devoted to the investigations of the concrete 
physical systems. Here we discuss how the general theory in the first half is used in the 
analysis of the concrete physical systems in the latter half. 

(A) In subsection 2.1, we gave the definitions of all the types of gauge invariant pertur- 
bation variables. The numerical data obtained by the cosmological observations are related 
with the gauge invariant perturbation variables, since both quantities do not depend on 
how we set up the spacetime coordinate system. Therefore it is desirable to express all the 
physical laws in the form closed by the gauge invariant quantities only. By the subsection 
2.1, it is guaranteed that the perturbation variables given by operating D{p) derivatives 
on the scalar like objects such as In a, sa '■= Pa/ P are gauge invariant. In the analysis of 
the concrete physical systems in the latter half of the present paper, the adiabatic per- 
turbation variable D(p)"lna, the entropic perturbation variables D{p)"^SA sue used as the 
independent variables and the total energy density p is used as the evolution parameter. 
This formalism by the A/E decomposition given in the subsection 3.2 and the p philosophy 
proposed in the subsection 3.3 is very useful for the physical interpretations of the results 
obtained in the latter half of the paper. 

(B) In the subsection 2.2, we construct the metric junction formalism as the method 
treating the sudden change of the equation of state. The metric junction formalism given 
in the subsection 2.2 is used in the analysis of the concrete physical systems in the latter 
half of the paper, since they contain the transitions such as the slow rolling oscillatory 
transition and the reheating transition. Since it is proven in the subsection 2.2 that our 
A/E perturbation variables D{p)"' In a, D{p)^SA are continuous across the matching surface 
defined by p = const, this set of variables must be useful in the research of the transitions 
even if we assume that the transition do not occur instantly. 

(C) Based on the p philosophy proposed in the subsection 3.3, in section 4 the evolutions 
of In a, Sa are described as the functions of the evolution parameter p. The evolutions of our 
A/E perturbation variables D(p)"'lna, D{p)"'SA are given by operating D{p) derivatives on 
the solutions of In a, sa expressed as the functions of p. Since D{p)f{p) = for an arbitrary 
function of p; f{p), the p dependences of In a, sa are directly reflected to the p dependences 
of D{p)"'\na, D{p)"'Sa- For this reason, by the formalism prepared in section 3, in section 
4 we can manifestly clarify the time evolutions in which the entropic perturbation of the 
energetically subdominant component makes the adiabatic perturbation D{p)"'lna grow. 
We can show that D{p)"'lna, the non-Gaussanity /tvl grow considerably and transiently 
when the energy ratio S2 of the energetically subdominant component which has the g 
factor smaller than the g factor of the dominant component pi and which is supported by 
the extremely small scalar field expectation value, begins to increase. Since until now as 
for D{p)"'lna only expression of the final state when the growth of D{p)"'lna has already 
ended, has been given, the p philosophy given in section 3 is superior in that the time 
evolution and the mechanism of the growth of the adiabatic perturbation variable called 
the Bardeen parameters D(p)"lna, the non-Gaussianity /tvl can be described. 

In the following paragraphs, the two points which was not treated in the first section 
are discussed. The first point is on the existence of more than two sources of cosmological 
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perturbations. The second point is on the modehng of the abrupt cosmic evolution by the 
metric junction across the spacehke hypersurface. 

(l)The present amphtude of the Bardeen parameter Cn(p) can be decomposed as Cn(p) = 
Cnsi(p) +Cnent(p)- Cnsi(p) IS the compoueut generated from the adiabatic mode in the slow 
rolling phase, and Cnent(p) is the adiabatic component generated from the entropy mode 
in the slow rolling phase by successive universe evolution. In many excellent papers, many 
authors insisted that a significant large nonlinearity can be generated in the inhomogeneous 
end of the inflation pHI] [291], in the modulated reheating [31]], in the curvaton scenario [32] 
] and in the vacuum dominated inflation [33] ] [M] ] . Unfortunately in the partial studies, 
without any plausible reasons it is assumed that Cnsi(p) is negligibly small compared with 
Cnent(p) and the non-Gaussianity parameter /tvl is calculated from Cnent(p) only. However 
the infiaton which drives the universe expansion enough to solve the horizon problem, 
the flatness problem, but do not generate any cosmo logical perturbations, does not exist. 
In this point of view, it is wonderful that the authors of the paper [32] ] tried to treat 
the contribution of the infiaton Cnsi(p) and the contribution of the curvaton Cnent(p) with 
equal importance from the standpoint of the mixed scenario. In this paper, we investigate 
whether a significant non-Gaussianity f^i is generated in the successive universe evolution 
taking into account the cosmological perturbations generated in the slow rolling phase. 
When we analyze the systems where more than two factors are concerned, the exponent 
evaluation method presented in this paper is very efficient. 

(2) In the early universe, there exists a period before which and after which the dynam- 
ical behaviors of each component are very different. As for the scalar field with mass m, 
while m -^ H, where H is the Hubble parameter, holds, the scalar field is in the slow 
rolling phase when its energy density changes mildly compared with the cosmic expansion, 
and while H <^m holds, the scalar field is in the oscillatory phase when its energy density 
behaves like a dust fiuid.[3|] [6|] [7]] [8]] In the reheating [H] [8|], the energy of the oscillatory 
scalar field behaving like a dust fiuid is transformed into that of the radiation fiuid. In 
the hybrid infiation, the energy of the slow rolling scalar field is transformed into that of 
the oscillatory scalar field and into that of the radiation fiuid on the bifurcation set. Such 
phase transitions are quite complicated and the completely rigorous mathematical treat- 
ment is beyond our scope. For example, we consider the slow rolling-oscillatory transition. 
In the m '^ H region and in the H '^ m region, the expansion schemes investigating 
the dynamical behaviors of the scalar field can be developed with m/H and H/m as the 
expansion parameter, respectively. P ] [3] [SI] [3] [3] But at the transition period H r^ m, 
any expansion schemes cannot be developed because of no expansion parameters. How- 
ever, in spite of complicated behaviors at the transition, the period of the transition can be 
thought to be short compared with the periods before and after the transition characterized 
hj m <^ H , H <^ m, respectively. Therefore we think the transitions as the instantaneous 
transient phenomena and may treat such transitions as the metric junctions across the 
spacelike hypersurfaces. In the above reasons, in our paper, the metric junction formula- 
tion on the cosmological perturbations in the long wavelength limit, linear and nonlinear, 
are constructed. On the spacelike hypersurface defined by H = m, the spacetime governed 
by the slow rolling scalar field and the spacetime governed by the oscillatory scalar field is 
connected. In case of reheating. On the spacelike hypersurface defined by if = F where F 
is the decay constant of the scalar field, the spacetime governed by the oscillatory scalar 



27 



field and the spacetime governed by the radiation fluid is connected. 
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§A Proof of Proposition 8 

Solving the matching condition C{t,x,X) = with respect to t gives t = t{x,X). By 
differentiating C{t, x, X) = with respect to A, x'', we obtain 

dt _ IdC 9t _ 1 

dX~~d^' d^'c' ^ ^ 

Differentiating [5*]^ = with respect to A, x^ gives 



dS ■dt^^ „ r. . ^dt 

dX^ dX 



ox' 



= 0, 
From the two sets of equations, we obtain [D{C)S]'t = 0, [Di{C)S]'t = 0. 



(A.2) 



§B the junction condition in the linear perturbation 
theory 

In this section, we consider the metric junction across the matching hypersurface charac- 
terized by C = where C is the scalar quantity in the linear perturbation theory, since the 
previous papers treating this problem [21]] [25]] [26]] sometimes contain some typographi- 
cal errors, derive the matching conditions without keeping the gauge invariance completely 
and are lacking in the physical interpretations from the viewpoint of the long wavelength 
limit. The notaions used in this section are based on the papers. [17]] [5]] [8]] The metric 
tensor is given by 

^00 = -{1 + 2AY), (B.l) 

go^ = -aBY,, (B.2) 

~g,, = a^[il + 2HLY)6,,+2HTY,j], (B.3) 

where Y, Yi and Yij are harmonic scalar, vector and tensor for a scalar perturbation with 
wavenumber k: 

where k"^ := ^^ fcj/cj. The energy momentum tensor of the total system is given by 

T^j^i, = {p + P)u^u^ + Pg^^, (B.5) 
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where p, P and u^ are the energy density, the pressure, and the four velocity of the total 
system. For the scalar quantities S = {p,P), S is expanded as S* = S* + SSY, and the four 
velocity u^ is written by 

uo = -{1 + AY), (B.6) 

Ui = a{v-B)Yi. (B.7) 

We define the gauge dependent geometrical quantities as 

n:=HL + ^Ht, ag := ^Ht - B. (B.8) 

For the scalar quantity S*, 

DS:=5S-^n (B.9) 

H 

is gauge invariant. For the four velocity -u^, the variable defined by 

Z:=n-^(v-B) (B.IO) 

k 



$:=7^-^-CT3 (B.ir 



is gauge invariant. The Newtonian potential $ defined by 

aH 

T 

is gauge invariant. 

We consider the metric junction across the hypersurface S defined by 

a;° = tx+5ZV, x' = y' + 6ZY\ (B.12) 

which connects the future spacetime Ai+ and the past spacetime A^_. From 6Z^, 6Z, we 
can define two gauge invariant quantities as 

0° := (5Z° + 47^, (B.13) 

(p := 5Z-\ht. (B.14) 

k 

We consider the case where the matching hypersurface S is defined by C = 0. In this case 

The normal vector of the matching hypersurface which points from A^_ to A^_|_ is given 
by 



n 



sgn{f)[-r''dpfd,f]-^'^d,f. (B.16) 



We define the intrinsic metric, the extrinsic curvature, the intrinsic energy momentum 
tensor by 



T 

f 



efe}'(5f^^ + n^n^), (B.17) 

gfg^^V^n,, (B.18) 

fi^fi%,, (B.19) 

n^e'^f^u, (B.20) 
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where ef := dx^ jdy^ . These quantities of the matching hypersurface E defined by (7 = 
can be written as 



9ii 



a' { 5,, + 25,,Y ( -H^ + ^'^ ) + '^M-k<P) \ , 



(B.21) 



a'HSi 



+ 



-2a''Hk^-e^-^^\Y, 



u 



C H 



ij: 



T 

f 



. f.DC k 



(B.22) 
(B.23) 
(B.24) 



where right hand sides are written in the gauge invariant form. The metric junction 
conditions across the matching hypersurface S are given by 

[g,,]+ = [k,,]t = [T„„]+ = [fn,]t = 0, (B.25) 

which yield the matching conditions in the long wavelength limit: in the background level, 

H+ = [H]t = [p]l = 0, (B.26) 

and in the perturbation level, 



0, 



DC 



[k'^^t = 0, 

-toe + Dp 
C 



^ ^ C H 



0. 



(B.27) 
(B.28) 

(B.29) 
(B.30) 



Owing to our previous paper |8]], in the long wavelength limit, the solution of DS where 
S is the scalar quantity is given by 



DS = DS^ + —C I dt^ 



(B.31) 



where 



DS^:-- 



dS S da 



dC^ a dC^ 



(B.32) 
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where dS/dC-,,, da/dCi, are the derivatives of the background quantities 5", a with respect 
to the solution constant C^ and c is a constant characterizing the adiabatic decaying mode. 
The solution of the Newtonian potential $ is given by 

e<^ = —c + 0{k^). (B.33) 

a 



Therefore (Ib:28|) . (E29]) give 



[c]+ = 0, (B.34) 

[D{C)at = [D{C)p]t = 0, (B.35) 



where as for scalar like object S 



^<^)^ - f - If- '^■^«) 



(IB.34|) . (IB.35|) are consistent with the metric junction conditions of the full nonlinear 



gradient expansion case represented by f l2.66p and the Proposition 8. 

§C the analyses of the evolution of the multiple vac- 
uum dominated scalar fields by the r function 

The T function was presented as the method of analyzing the evolution of the multiple 
scalar fields. P] Under the slow rolling approximation, the evolution of the scalar fields is 
described by 

d(pa L^ (r..-. 

dN K'^Ud(t)a ■ 

which are decomposed as 

* ^" '"' .^U, (C.2) 



dr d(f)a dr 

introducing the r function as the new evolution parameter. P ] It is much easier to treat 
the new evolution equations (IC.2I) than to treat the original evolution equations ( JC.lll for 
many cases. 

We consider the vacuum dominated case given by 

U = Uo + J2 lva<l>l. (C.3) 

a 

In this case, the evolution of (pa is given by 

0a = 0,(O)exp(-r7,r). (C.4) 

By using the r function as the evolution parameter, the Bardeen parameter Cn(p) are given 

by 
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where the subscripts A, r are interpreted as the derivatives with respect to A, r, respectively; 
for example 

U. := llu. (C.6) 

Concretely Ci{p), C2(p) are given by 

-^Ci(p) = [drU.-^U, (C.7) 

^C2(P) = fdTUxx-^ + ^(-UxxUr + 2UxUxr-^Urry (C.8) 

and the coefficients A^'^, Nat are given by the expansions 

Up) = 1^N^^^, Up) = I^N^.^^ ^, (C.9) 

a ab 

where we assume that all the nonlinear perturbations at the ffist horizon crossing in the 
inflationary expansion are vanishing: (i"0a(O)/(iA" = (n > 2). By using A{n^k) defined 

by 

A{n,k) := 5^r/>2(o)exp(-2fc77,r), (CIO) 



and collecting the leading order terms with respect to Uq, we obtain 

(c.ii) 



. 1 

JNL 



K^o 



A{2,1)A{3,3) ^ A{3,2) ^A{3,1) 



A(2,2)2 A(2,2) A(2,l) 



We consider the system where two scalar fields evolve. Since only the first term in (IC.lip 
can change the exponent for the moving r, we concentrate on the first term written by 
(Inl)! from now on. {fNLli can be written by 

where 

Vr-=—, 0r := -^, e^ :=exp{-2(r/2-^i)r}. (C.13) 

Vi 01 (0) 

In order that the infiation can solve the horizon problem, we assume rja/n'^UQ ~ 10~^. 
First we consider the case ?7i > 0, 772 < 0. we put 

0i(O) = llO™, 02(0) = -10-', e, = 10^ (C.14) 

K K 



then we get 



{fNL)l - iO (1 + 10~2a+m)+2p)2 " y'^'^^) 
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The exponent evaluation method gives 
P <Pi 

Pl <P <P2 
P2 <P <P3 
P3<P 



where 



2 

Pi-=^{i + m), p2 



{fNL)i =lo-^ 


(C.16) 


ifNLh = - lo3^-^('+")~^ 


(C.17) 


Unl)! = - lo-^+^('+")-^ 


(C.18) 


(/7VL)l=-10-^ 


(C.19) 


m), ps:=2{l + m). 


(C.20) 



The above evaluations agree completely with fl4.59p . fl4.60p . fl4.6ip by taking the correspon- 
dence p = k + 2m — 2. 

Next we consider the case where rja < 0, rji > rj2- We put 

0i(O) = -10^'", 02(O) = -1O~', er = W, (C.21) 

K K 

where m < I, then we obtain 

[fNLh - -lU ^^ ^ ^Q_2{l-m)+2py • ^'^•^^^ 

The exponent evaluation method gives 



<PI 








ifNL)i = - lo-^ 


(C.23) 


<P <P2 








{fNLh = - lo3^-2('-")-^ 


(C.24) 


<P<P3 








Unl)! = - lo-^+^('-")-^ 


(C.25) 


<p 








(/7VL)l=-10-^ 


(C.26) 


Pi--=b- 


-m), 


P2 ■■ 


= {1- 


-m), p3:=2(/-m). 


(C.27) 



where 



The above evaluations agree completely with fl4.69p . (14 .70 p . fl4.7ip by taking the correspon- 
dence p = k — 2m — 2. 

References 

[1] Bardeen, J.M., Phys. Rev. D 22, 1882 (1980). 

[2] Gordon, C, Wands, D., Bassett, B.A. and Maartens, R., Phys. Rev D63 023506 
(2000). 

[3] Kodama, H. and Hamazaki, T., Prog. Theor. Phys. 96,949-970 (1996). 

[4] Hamazaki, T. and Kodama, H., Prog. Theor. Phys. 96,1123-1146 (1996). 

[5] Kodama, H. and Hamazaki, T., Phys. Rev. D57, 7177-7185 (1998). 

[6] Hamazaki, T., Phys. Rev. D 66, 023529 (2002). 

33 



[7] Hamazaki, T., Nucl. Phys. B 698,335-385 (2004). 

[8] Hamazaki, T., Nucl. Phys. B 791,20-59 (2008). 

[9] Hamazaki, T., Phys. Rev. D 78, 103513 (2008). 

[10] Taruya, A. and Nambu, Y., Phys. Lett. B428 37-43 (1998). 

[11] Sasaki, M. and Tanaka, T , Prog. Theor. Phys. 99,763-782 (1998). 

[12] Wands, D., Mahk, K.A., Lyth, D.H. and Liddle, A.R. Phys. Rev. D62, 043527 (2000). 

[13] Rigopoulos, G.I. and Shellard, E.P.S. Phys. Rev. D68, 123518 (2003). 

[14] Lyth, D.H., Mahk, K.A. and Sasaki, M. JCAP. 0505, 004 (2005). 

[15] Nakamura, K. Phys. Rev. D74, 101301 (2006). 

[16] Tanaka, Y. and Sasaki, M. , Prog. Theor. Phys. 117,633 (2007). 

[17] Kodama, H. and Sasaki, M., Prog. Theor. Phys. Suppl. 78, 1-166 (1984). 

[18] Kodama, H. and Sasaki, M., Int. J. Mod. Phys. A2, 491 (1987). 

[19] Mukhanov, V.F., Feldman, H.A. and Brandenberger, R.H. Phys. Rep. 215, 203 (1992). 

[20] Langlois, D. , Vernizzi, F. and Wands, D., JCAP0812:004 (2008) 

[21] Takahashi, T., Yamaguchi, M., Yokoyama, J. and Yokoyama, S. Phys. Lett. B678:15 
(2009). 

[22] Takamizu, Y., Mukohyama, S., Sasaki, M. and Tanaka, Y. JCAP1006:019 (2010). 

[23] Israel, W. Nuovo. Cim. 44B, 1 (1966). 

[24] Deruelle, N. and Mukhanov, V.F. Phys. Rev. D52, 5549 (1995). 

[25] Martin, J. and Schwarz, D.J. Phys. Rev. D57, 3302 (1998). 

[26] Copeland, E.J. and Wands, D. JCAP0706:014 (2007). 

[27] Mukohyama, S. Class. Quant. Grav. 17:4777 (2000). 

[28] Lyth, D.H. JCAP0511:006 (2005). 

[29] Sasaki, M. , Prog. Theor. Phys. 120,159 (2008). 

[30] Sasaki, Y. and Yokoyama, J. Phys. Rev. D44, 970 (1991). 

[31] Dvah, G., Gruzinov, A. and Zaldarriga, M., Phys. Rev. D 69, 023505 (2004). 

[32] Ichikawa, K., Suyama, T, Takahashi, T and Yamaguchi, M. Phys. Rev. D78, 023513 
(2008). 

[33] Alabidi, L. JCAP0610:015 (2006). 

[34] Byrnes, C.T., Choi, K.Y. and Hall, L.M.H JCAP0902:017 (2009). 

34 



[35] Komatsu, E. and Spergel, D.N. Phys. Rev. D63, 063002 (2001). 

[36] Malik, K.A. and Wands, D., Class. Quant. Grav. 21, L65 (2004). 

[37] Malik, K.A., JCAP. 0511, 005 (2005). 

[38] Shibata, M. and Sasaki, M. Phys. Rev. D60, 084002 (1999). 



35 



